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A search for — >■ fi'^ fJ.~ and — >■ fJ-^ fi~ decays is performed using 1.0 fb~^ of pp collision data 
collected at = 7 TeV with the LHCb experiment at the Large Hadron Collider. For both decays 
the number of observed events is consistent with expectation from background and Standard Model 
signal predictions. Upper limits on the branching fractions are determined to be B{B^^ — > fj,'^fj,~) 
< 4.5(3.8) X 10"'' and B{B° At'^M") < 1-0(0.81) x 10"^ at 95% (90%) confidence level. 



Flavor changing neutral current (FCNC) processes are 
highly suppressed in the Standard Model (SM) and thus 
constitute a stringent test of the current description of 
particle physics. Precise predictions of the branching 
fractions of the FCNC decays 5° M^/^^ and B° 
fi+^i-, BiB° ^ pt+fL') = (3.2±0.2) X 10-9 and B{B^ 
11+ H ) = (0.10 ± 0.01) X 10-9 [1113 make these modes 
powerful probes in the search for deviations from the 
SM, as contributions from new processes or new heavy 
particles can significantly modify these values. Previous 
searches [3HB] already constrain possible deviations from 
the SM predictions, with the lowest published limits from 
the LHCb collaboration: S(B° A^^M") < 1-4 x 10"* 
and B{B° fi+H^) < 3.2 x 10"^ at 95% Confidence 
Level (CL). 

In this Letter, we report an analysis of the pp collision 
data recorded in 2011 by the LHCb experiment corre- 
sponding to an integrated luminosity of 1.0 fb-^. This 
dataset includes the 0.37 fb-^ used in the previous anal- 
ysis f^. In addition to the larger dataset, improvements 
include an updated event selection, an optimized bin- 
ning in the discriminating variables, and a reduction of 
the peaking background. The data already analyzed in 
Rcf. [6' were reprocessed and, to avoid any potential bias. 



all the events in the signal region were blinded until all 
the analysis choices were finalized. 

The LHCb detector [7] is a single-arm forward spec- 
trometer covering the pseudo-rapidity range 2 < 77 < 5. 
The detector includes a high precision tracking system 
consisting of a silicon-strip vertex detector, a large- 
area silicon-strip detector located upstream of a dipole 
magnet with a bending power of about 4Tm, and 
three stations of silicon-strip detectors and straw drift 
tubes placed downstream. The combined tracking sys- 
tem has a momentum resolution A.p/p that varies from 
0.4% at 5GeV/c to 0.6% at lOOGeV/c. Two ring- 
imaging Cherenkov detectors (RICH) are used to identify 
charged particles. Photon, electron and hadron candi- 
dates are identified by a calorimeter system consisting 
of scintillating-pad and pre-shower detectors, an electro- 
magnetic calorimeter and a hadronic calorimeter. Muons 
are identified by alternating layers of iron and niultiwire 
proportional chambers. 

The trigger consists of a hardware stage, based on in- 
formation from the calorimeter and muon systems, fol- 
lowed by a software stage (HLT) that applies a full 
event reconstruction. Events with muon final states 
are triggered using two hardware trigger decisions: the 
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single-muon decision (one muon candidate with trans- 
verse momentum > 1.5GeV/c), and the di-muon deci- 
sion (two muon candidates with pt,i and pT,2 such that 
y/PTXPr^ > 1.3GeV/c). All tracks in the HLT are re- 
quired to have a, pt > 0.5 GeV/c. The single muon trig- 
ger decision in the HLT selects tracks with an impact 
parameter IP > 0.1mm and px > l.OGeV/c. The di- 
muon trigger decision requires fJ.^ fJ-~ pairs with an in- 
variant mass m^fj_ > 4700 MeV/c^. Another trigger de- 
cision, designed to select J/ip mesons, requires 2970 < 
m^^ < 3210MeV/c^. Events with purely hadronic final 
states are triggered by the hardware trigger if there is a 
calorimeter cluster with transverse energy Et > 3.5 GeV. 
HLT trigger decisions selecting generic 6-hadrons decays 
provide high efficiency for such final states. 

The S"^^ fJ.^ii~ selection requires two high qual- 
ity muon candidates displaced with respect to any pri- 
mary pp interation point (primary vertex, PV). The di- 
muon secondary vertex (SV) is required to be well mea- 
sured (with a per degree of freedom smaller than 9.0), 
downstream, and separated from the PV by a distance- 
of- flight significance greater than 15. When more than 
one PV is reconstructed, the one giving the minimum IP 
significance for the B candidate is chosen. Only candi- 
dates with IP/(t(IP) < 5 are kept. Combinations with 
poorly reconstructed tracks are removed by requiring 
p < 500GeV/c and 0.25 < pt < 40GeV/c for all tracks 
from the selected candidates. Only B candidates with 
decay times smaller than 9 x t{B^) [S] are kept. Fi- 
nally, according to the simulation, approximately 90% of 
di-muon candidates coming from elastic di-photon pro- 
duction are removed by requiring a minimum of the 
B candidate of 500MeV/c. The surviving background 
mainly comprises random combinations of muons from 
semileptonic b-hadrons decays (bb — > X, where X 

is any other set of particles) . 

Three channels, B+ J/ipK"^, B° and 
B^ X+TT" (inclusion of charged conjugated processes 
is implied throughout this Letter) serve as normalization 
modes. The first two have trigger and muon identifi- 
cation efhciencies similar to those of the signal, but a 
different number of tracks in the final state. The third 
channel has a similar topology, but is selected by differ- 
ent triggers. The selection of these channels is designed 
to be as similar as possible to that of the signal to reduce 
the impact of common systematic uncertainties. An in- 
clusive B'^^s^ — > h^h'~ sample (where h,h' can be a pion 
or a kaon) is the main control sample. The selection 
is the same as for -> fi^p^ signal candidates, ex- 
cept for the muon identification requirement. To ensure 
similar selection efficiencies for the B° K^tt" and 
^{s) ~^ channels, tracks from the 5° K^n^ 

decay are required to be in the muon detector accep- 
tance. The J/ip — >■ decay in the i?+ — ^ J/ipK'^ 
and — J/^^P normalization channels is also selected 



as — >■ ^^^^ signal, except for the requirements on 
its IP and mass. Kaon candidates are required to be 
identified by the RICH detectors and to pass IP selec- 
tion criteria. 

A multivariate selection (MVS), based on a boosted de- 
cision tree removes 80% of the residual background, 
while retaining 92 % of the signal. Applying this selec- 
tion improves the performance of the main multivariate 
algorithm described below. The six variables entering 
the MVS, ordered by their background rejection power, 
are: the angle between the direction of the momentum of 
the B candidate and the direction defined by the vector 
joining the secondary and the primary vertices, the B 
candidate IP and its vertex x^, the minimum IP of the 
muons with respect to any PV, the minimum distance 
between the two daughter tracks and the of the SV. 
The B^^-j — ?► h^h'^ mass sidebands have been used to 
check that the distribution of the MVS output is similar 
for data and simulation. The same selection is applied 
(using, when necessary, slightly modified variable defi- 
nitions) to the normalization samples. The efficiencies 
for the signal and the normalization samples are equal 
within 0.2 % according to the simulation. 

In total, 17321 muon pairs with invariant mass be- 
tween 4900 and 6000MeV/c^ pass the trigger and se- 
lection requirements. Given the measured bb cross- 
section [To] and assuming SM rates, this data sam- 
ple is expected to contain 11.6 5° — > /i^/i^ and 1.3 
B° — > decays. 

The selected candidates are classified in a binned two- 
dimensional space formed by the di-muon invariant mass 
and the output of another boosted decision tree, de- 
scribed in detail below. In the following we employ 
"BDT" to indicate the algorithm or its ouput, depending 
on the context. 

The invariant mass line shape of the signal events is 
described by a Crystal Ball function |TT]. The peak val- 
ues for the B^ and B^ mesons, m^o and mso, are ob- 
tained from the 5° — ?> K^K^ and B'^ — > K^tt^ sam- 
ples [12]. The resolutions are extracted from data with 
a power-law interpolation between the measured resolu- 
tions of charmonium and bottomonium resonances de- 
caying into two muons. Each resonance is fitted with 
the sum of two Crystal Ball functions with common 
mean values and resolutions, but different parameters 
describing the tails. The results of the interpolation 
at TO^o and niBO are ct(tobo) — 24.8 ± 0.8MeV/c^ and 
ct(tobo) = 24.3±0.7MeV/c^. They are in agreement with 
those found using B^ — >■ K~^Tr~ and 5° — >■ K~^K~ exclu- 
sive decays. The transition point of the radiative tail is 
obtained from simulated B'^ events re-weighted 

to reproduce the mass resolution measured in data. 

Geometrical and kinematic information not fully ex- 
ploited in the selection is combined via the BDT for which 
nine variables are employed [6j. Ordered by their back- 
ground rejection power, they are: the B candidate IP, the 
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minimum IP significance, the sum of the degrees of isola- 
tion of the muons (the number of good two-track vertices 
a muon can make with other tracks in the event), the B 
candidate decay time, p^, and degree of isolation [13], 
the distance of closest approach between the two muons, 
the minimum px of the muons, and the cosine of the an- 
gle between the muon momentum in the di-muon rest 
frame and the vector perpendicular to the B candidate 
momentum and to the beam axis. No data were used for 
the choice of the variables and the subsequent training of 
the BDT, to avoid biasing the results. Instead the BDT 
was trained using simulated samples {B'^^^ — > fJ.^n'' for 
signal and bb — for background). The BDT out- 
put is independent of the invariant mass for signal inside 
the search window. It is defined such that for the signal 
it is approximately uniformly distributed between zero 
and one, while for the background it peaks at zero. 

The probability for a signal event to have a given 
BDT value is obtained from data using an inclusive 



BDT and invariant mass bin inside the signal regions is 
determined from data by fitting to an exponential func- 
tion events in the mass sidebands defined by [4900, 5000] 
MeV/c^ and [mso +60 MeV/c^ 6000 MeV/c^ ]. The 
systematic uncertainty on the estimated number of com- 
binatorial background events is computed by fluctuating 
with a Poissonian distribution the number of events mea- 
sured in the sidebands, and by varying within ±la the 
value of the exponent. As a cross-check, another model, 
the sum of two exponential functions, has been used to 
fit the events in different ranges of sidebands providing 
consistent background estimates inside the signal regions. 
An additional systematic uncertainty is introduced where 
the yields in the signal regions differ by more than 1 a be- 
tween the fit models. 

Peaking backgrounds from S"^^ h^h'~ events have 
been evaluated by folding the K fi and tt — > /i 
misidentification rates extracted from a 1?° — > A'^7r+ 
sample from data in bins of p and pt into the spec- 



5° ^ h+h'- sample. Only events triggered indepen- trum of selected simulated B^^ h+h' events. The 



dently of the presence of any track from the signal candi- 
dates are considered. The number of B^^^ — >■ h^h'" sig- 
nal events in each BDT bin is determined by fitting the 
hh' invariant mass distribution. The maximum spread 
in the fractions of the yields going into each bin, ob- 
tained by fitting the same dataset with different signal 
and background models, is used to evaluate the system- 
atic uncertainty on the signal BDT probability distribu- 
tion function _6J. 

The binning of the BDT and invariant mass distri- 
butions is re-optimized with respect to Ref. iSj, using 
simulation, to maximize the separation between the me- 
dian of the test statistic distribution expected for back- 
ground and SM B^ — signal, and that expected 
for background only. The chosen number and size of the 
bins are a compromise between maximizing the number 
of bins and the necessity to have enough — > h^h'" 
events to calibrate the B^ — )• BDT and enough 

background in the mass sidebands (see below) in each 
bin to estimate the combinatorial background in the 
S° and i?° mass regions. The BDT range is thus di- 
vided into eight bins (see Table [l]) and the invariant 
mass range into nine bins with boundaries are defined by 
m^o ± 18, 30, 36, 48, 60 MeV/c^. This binning improves 
the test statistic separation by about 14% at the SM 
rate with respect to Ref. [6 ; over 97 % of this separation 
comes from the bins with BDT> 0.5. 

We select events in the invariant mass range 
4900 MeV/c^, 6000 MeV/c^]. The boundaries of the signal 
regions are defined as m^o ^ ± 60 MeV/ . The low-mass 
sideband is potentially polluted by cascading b — >■ c/ii/ — ?> 
HfiX decays below 4900 MeV/ (? and peaking background 
from — > h^h' candidates with the two hadrons 
misidentified as muons above 5000 MeV/ c^. The num- 
ber of expected combinatorial background events in each 



mass line shape of the peaking background is ob- 
tained from a simulated sample of doubly-misidentified 
h+h'^ events. In total, O.btal (2.6^^:4) doubly- 
h^h' events are expected in the 
B^ {B^) signal mass windows. The contributions of 
5+ — ^ J/tp{fi'^ and B^ — >• exclusive de- 

cays have been found to be negligible with respect to the 
combinatorial and i?"^^ — > h^h' backgrounds. 

The 5° — 7> f^^lJ-^ and B'^ — > ji^ yields are translated 
into branching fractions using 



misidentified B"^^ 



/n 



d{s 



An 



(1) 



where fd(s) and /norm are the probabilities that a b quark 
fragments into a B'^^^ and into the hadron involved in the 
given normalization mode respectively. We use fs/ fd = 
0.267t!]:[!^J \XM and we assume fd = /„• With Bnorm 
we indicate the branching fraction and with A'norm the 
number of signal events in the normalization channel ob- 
tained from a fit to the invariant mass distribution. The 
efficiency esig(norm) for the signal (normalization channel) 
is the product of the reconstruction efficiency of all the 
final state particles of the decay including the geomet- 
ric acceptance of the detector, the selection efficiency for 
reconstructed events, and the trigger efficiency for recon- 
structed and selected events. The ratio of acceptance and 
reconstruction efficiencies are computed using the Monte 
Carlo simulation. The differences between the simulation 
and data are included as systematic uncertainties. The 
selection efficiencies are determined using Monte Carlo 
simulation and cross-checked with data. Reweighting 
techniques have been used for all the Monte Carlo dis- 
tributions that do not match those from data. The trig- 
ger efficiency is evaluated with data driven techniques. 
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Finally, A^^o is the number of observed signal 

events. The observed numbers of — >■ J /ipK'^ , — >■ 
J/ij}(t) and K+n- candidates are 340 100 ± 4500, 

19 040 ± 160 and 10 120 ± 920, respectively. The three 
normalization factors are in agreement within the uncer- 
tainties and their weighted average, taking correlations 
into account, gives = (3.19 ± 0.28) x 10"^° 

and a^°o™^+^_ = (8.38 ± 0.39) x 10~". 

For each bin in the two-dimensional space formed by 
the invariant mass and the BDT we count the number 
of candidates observed in the data, and compute the ex- 
pected number of signal and background events. 

The systematic uncertainties in the background and 
signal predictions in each bin are computed by fluctu- 
ating the mass and BDT shapes and the normalization 
factors along the Gaussian distributions defined by their 
associated uncertainties. The inclusion of the systematic 
uncertainties increases the — )■ /x^/Lt" and — )■ ^"""/i^ 
upper limits by less than ^ 5%. 

The results for 5° — >• and B^ — > decays, 

integrated over all mass bins in the corresponding signal 
region, are summarized in Table [l] The distribution of 
the invariant mass for BDT>0.5 is shown in Fig. [l] for 
B° — >■ fi^n^ and i?" — > candidates. 




qI I & 111*1 I 

5350 5400 5250 5300 

m,,(MeV/c') m,,(MeV/c') 



FIG. 1. Distribution of selected candidates (black points) 
in the (left) — >■ /i^/i^ and (right) M^M^ mass 

window for BDT>0.5, and expectations for, from the top, 
^(s) SM signal (gray), combinatorial background 

(light gray), B^^j — > K^h'~ background (black), and cross- 
feed of the two modes (dark gray). The hatched area depicts 
the uncertainty on the sum of the expected contributions. 

The compatibility of the observed distribution of 
events with that expected for a given branching frac- 
tion hypothesis is computed using the CLg method [IS] . 
The method provides CLg+b, a measure of the com- 
patibility of the observed distribution with the signal 
plus background hypothesis, CLb, a measure of the 
compatibility with the background-only hypothesis, and 
CLg = CLs+b/CLb. 



The expected and observed CL^ values are shown in 
Fig. [2] for the i?^ — >■ /x+^i" and 5° — >■ channels, 
each as a function of the assumed branching fraction. 
The expected and measured limits for i?^ — )■ /i^/i^ and 
B° ^i+tjL- at 90 % and 95 % CL are shown in Table |ll] 
The expected limits are computed allowing the presence 
of -B^^j — ?► /x^M~ events according to the SM branching 
fractions, including cross-feed between the two modes. 

The comparison of the distributions of observed 
events and expected background events results in a p- 
value (1-CLb) of 18% (60%) for the ^ 
[B^ — > /i+/i^) decay, where the CLb values are those cor- 
responding to CLg+b = 0.5. 

A simultaneous unbinned likelihood fit to the mass pro- 
jections in the eight BDT bins has been performed to 
determine the B^ — )• jJL^ n,~ branching fraction. The sig- 
nal fractional yields in BDT bins are constrained to the 
BDT fractions calibrated with the B^ — > h^h'^ sam- 

pie. The fit gives B{B° -> m+m") = (0.8ti:3) x 10"^ 
where the central value is extracted from the maximum 
of the logarithm of the profile likelihood and the uncer- 
tainty reflects the interval corresponding to a change of 
0.5. Taking the result of the fit as a posterior, with a 
positive branching fraction as a flat prior, the probabil- 
ity for a measured value to fall between zero and the SM 
expectation is 82 %, according to the simulation. The 
one-sided 90%, 95% CL hmits, and the compatibility 
with the SM predictions obtained from the likelihood, are 
in agreement with the CLg results. The results of a fully 
unbinned likelihood fit method are in agreement within 
uncorrelated systematic uncertainties. The largest sys- 
tematic uncertainty is due to the parametrization of the 
combinatorial background BDT. 

In summary, a search for the rare decays — > /i^/i^ 
and B'^ — > /i^/i^ has been performed on a data sam- 
ple corresponding to an integrated luminosity of l.Ofb"^. 
These results supersede those of our previous publica- 
tion [6] and are statistically independent of those ob- 
tained from data collected in 2010 [H]. The data are 
consistent with both the background-only hypothesis and 
the combined background plus SM signal expectation at 
the 1 a level. For these modes we set the most stringent 
upper limits to date: B{B^ — > fi'^f-i^) < 4.5 x 10^^ and 
B{B" < 1.03 X 10-9 at 95% CL. 
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FOM and NWO (The Netherlands); SCSR (Poland); 
ANCS (Romania) ; MinES of Russia and Rosatom (Rus- 
sia); MICINN, XuntaGal and GENCAT (Spain); SNSF 
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TABLE I. Expected combinatorial background, Bj'^j — > h^h' background, cross-feed, and signal events assuming SM pre- 
dictions, together with the number of observed events in the — > fJ.^ fJ-^ and B" — >■ fi"^ ii~ mass signal regions, in bins of 
BDT. 



Mode 



BDT bin 



0.0 - 0.25 0.25 - 0.4 0.4 - 0.5 0.5 - 0.6 0.6 - 0.7 0.7 - 0.1 



0.8-0.9 0.9-1.0 



^ jj. Exp. comb, bkg 



1889t^» 



57 



-11 
-11 



4.3 



-1-1.0 
1.0 



Exp. peak, bkg 0.124lg:°^^ OMStomt 0.0491^12 0.0451;;; 
Exp. signal 2.551;;;™ l-22t°;?° 0.971;; 
Observed 1818 39 12 



016 
012 



"•»D^-0.088 



q qr| + 0.92 
O.OU_Q gg 

0.0501^1^3 

i-ool";;^ 



1 06+°-^^ 

i.>J'J_0.46 

1 034+0109 



1 97+0.53 
-0.52 

0.0491°;°}^ 
1.18l°;i? 



44+0-'*! 

1 90+0.21 

^■^•^-0.21 



Exp. comb, bkg 
Exp. peak, bkg 
Exp. cross-feed 
Exp. signal 
Observed 



2OO3I43 
0.7ll»;re 

0.401°;}^ 
S00+0 0**6 

1904 



61 



1-12 
-11 



16.611? 



+0.71 
0.50 



4.711J 3.52lJ;J? 1.11 
0.355l°;Sf3 0.279tl:lll 0.249tt°oll 0.280l°;r2 0.264l03;io? 

0.136l°;»}I 0.158lg;°l? 0.164l0;0}? 



1 62+0™ 

n 97c;+0.108 
U.z I <J_o.060 



n iqQ+0.033 

n i4c;+0 027 
50 



n 1 KQ+0.023 
U.100_Q Q21 



0.115 
20 



■0.020 
0.017 



102+0014 

"•1"^-0.013 



iiq+0.017 



n 1 90+0.016 

U.lZO_Q^Q]_g 



0.187 
0.140; 



h0.022 
-0.020 

hO.019 
-0.017 



S4+0-53 

267+0 106 

104+0 036 
U.iy4_Q Q33 

145+0-030 
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B(B^^H- 

2. CLs as a function of the assumed B for (left) B'^ 



1.2 1.4 1.6 1.8 

,0 



0.2 0.4 0.6 0.8 1 

B(B"^H*H-) [10""] 

li~ and (right) — > fJ-^ IJ-^ decays. The long dashed black 
fj,'^lj,~ , if background and SM signal were observed, and for 



FIG. 2. UIjs as a tunction ot tne assumed O tor (lettj — >■ 
curves are the medians of the expected CLs distributions for B^ 

B^ — ^ fj,'^ , if background only was observed. The yellow areas cover, for each B, 34% of the expected CLs distribution on 
each side of its median. The solid blue curves are the observed CLs. The upper limits at 90% (95%) CL are indicated by the 
dotted (solid) horizontal lines in red (dark gray) for the observation and in gray for the expectation. 



TABLE II. Expected and observed limits on the ijO^j — >■ fi'^ /i 
branching fractions. 



ration, the CMS collaboration released the results of an 
updated search for these channels [IB] . 



Mode 




Limit 


at 90 % CL 


at 95 % CL 






Exp. bkg-hSM 


6.3 X 10"'' 


7.2 X 10"" 






Exp. bkg 


2.8 X 10"^ 


3.4 X 10^9 






Observed 


3.8 X 10-9 


4.5 X 10-9 


fiO^ 




Exp. bkg 


0.91 X 10"'^ 


1.1 X 10-^ 






Observed 


0.81 X lO"'* 


1.0 X lO-'' 
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